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ABSTRACT 

We present a comprehensive catalog of high signal-to-noise spectra obtained 
with the DEIMOS spectrograph on the Keck II telescope for a sample of 
F850LP<22.43 (AB) held spheroidal (E+SOs; 163) and bulge dominated disk 
(61) galaxies in the redshift range 0.2 < z < 1.2, selected on the basis of vi¬ 
sual morphology from the northern held of the Great Observatories Origins Deep 
Survey (GOODS-N). We discuss sample selection, photometric properties and 
spectral reduction. We derive scale-lengths, surface brightnesses and photometric 
inhomogeneities from the ACS data, and redshifts, stellar velocity dispersions, 

[O II] and HJ equivalent widths from the Keck spectroscopy. Using the pub¬ 
lished 2Ms Chandra X-ray catalog we identify active galactic nuclei to clarify 
the origin of emission lines seen in the Keck spectra. Only 2/13 [O II] emitting 
early-type galaxies are identihed as secure AGN on the basis of their X-ray emis¬ 
sion. Contrary to earlier suggestions, we hnd that most spheroidals containing 
‘blue cores’ are not associated with non-thermal nuclear activity. We examine 
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the zero point, tilt and scatter of the Fundamental Plane (FP) as a function 
of redshift and morphological properties, carefully accounting for luminosity- 
dependent biases via Montecarlo simulations. The evolution of the overall FP 
can be represented by a mean change in effective mass-to-light ratio given by 
< d\og{M/L- q)/ dz >= —0.72]'^Q;g5 ± 0.04. However, this evolution depends sig- 
nihcantly on the dynamical mass, being slower for larger masses as reported in a 
previous letter. In addition, we separately show the intrinsic scatter of the FP in¬ 
creases with redshift as d(rms(M/LB))/dz=0.040±0.015. Although these trends 
are consistent with single burst populations which formed aX Zf > 2 for high 
mass spheroidals and Zf ~ 1.2 for lower mass systems, a more realistic picture 
is that most of the stellar mass formed in all systems aX z > 2 with subsequent 
activity continuing to lower redshifts {z < 1.2). The fraction of stellar mass 
formed at recent times depend strongly on galactic mass, ranging from < 1% 
for masses above IO^^'^Mq to 20-40% below 10^^ Mq. Independent support for 
recent activity is provided by spectroscopic ([O II] emission, H(5) and photomet¬ 
ric (blue cores and broad-band colors) diagnostics. Via the analysis of a large 
sample with many independent diagnostics, we are able to reconcile previously 
disparate interpretations of the assembly history of held spheroidals. We discuss 
the implications of this measurement for the determination of the evolution of 
the number density of E-|-S0s galaxies, suggesting number density evolution of 
the morphologically selected population has occurred since 2 : ~1.2. 

Subject headings: cosmology: observations, galaxies: formation, galaxies: evolu¬ 
tion, morphologies 


1. Introduction 

The assembly history of held spheroidal galaxies (dehned here to include both ellipti¬ 
cals and lenticulars, in short E-|-S0s) remains an important issue in galaxy formation and 
particularly in the interface between theory and observation (see Treu 2004, McCarthy 2004 
for recent reviews). The overall assembly rate of stellar mass has been constrained both by 
measures of the cosmic star formation history (Madau et al. 1996, Blain et al. 1999, Steidel 
et al. 1999, Bunker et al. 2004) and, directly, via photometric estimates of the stellar mass 
density as a function of redshift (Brinchmann & Ellis 2000, Cohen 2001, Papovich et al. 
2001, Rudnick et al. 2003; Fontana et al. 2004). However, such global constraints integrate 
over important details relating to the formation history of individual morphological types. 

A key prediction of current semi-analytical models based on the cold dark matter 
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paradigm (e.g. Kauffmann 1996, Baugh et al. 1998, Somerville & Primack 1999) is the 
late assembly of massive spheroids from smaller fragments with signihcant consequences for 
the abundance and mass distribution of systems at ~l-2. Although the actual rate of 
assembly is sensitive to details of the algorithms used in the semi-analytical codes (see com¬ 
parisons in Benson et ah 2003), the observational data has provided only mixed support for 
the evolution expected. 

The simplest evolutionary test is based on measuring the density of spheroids at various 
redshifts. Brinchmann & Ellis (2000), and more recently Bundy, Ellis & Conselice (2005), 
found only a modest decline in the stellar mass density to z ~1 using morphologically-dehned 
samples (see also Im et ah 2002, Cross et al. 2004). These results appear to contrast with 
the claim by Bell et al. (2003), based on color-selected samples, that the luminosity density 
associated with red galaxies is fairly constant to ~1 suggesting a marked drop in the 
volume density of such objects. The utility of these tests depends sensitively on corrections 
for incompleteness and cosmic variance. 

An independent approach involves estimating the ongoing mass accretion rate for spheroidals, 
either via the presence of photometric abnormalities such as ‘blue cores’ (Menanteau et al. 
2001a,b, 2004, 2005), spectroscopic signatures of recent activity (Trager et al. 2000; Treu 
et al. 2002; van Dokkum & Ellis 2003) or deviations from passive trends predicted in the 
evolution of the Fundamental Plane (Djorgovski & Davis 1987, Dressier et al. 1987; hereafter 
FP; Treu et al. 1999, 2001b, 2002; van Dokkum & Ellis 2003; Gebhardt et al. 2003; van de 
Wei et al. 2004; Treu & Koopmans 2004; see also Kochanek et al. 2000; Rusin et al. 2003a; 
van de Ven et al. 2003; Rusin & Kochanek 2005). The degree to which the population is 
homogeneous in its passive evolution provides a very sensitive measure of the rate at which 
spheroidals may still be forming. Both tests are valuable but progress on the latter has, until 
now, been hindered by the small size of the spectroscopic samples available. 

This paper is motivated by extending the early work begun in this area via a comprehen¬ 
sive analysis of high quality Keck spectroscopy for 224 morphologically-selected spheroidals 
and bulge-dominated galaxies in the northern held of the Great Observatories Origins Deep 
Survey (GOODS). Associated with the Advance Gamera for Surveys (AGS) and spectroscopic 
data discussed here, is iP-band photometry from which stellar masses and optical-infrared 
colors have also been derived elsewhere (Bundy, Ellis & Gonselice, 2005). 

An overview of one of the main results of this observational program has already been 
presented in Treu et al. (2005; hereafter T05), where we discuss the evolution of dynamical 
mass-to-light ratios as a function of mass and environmental density in the context of hierar¬ 
chical models of galaxy formation. We hnd that mass-to-light ratios - and hence luminosity 
weighted stellar ages - depend sensitively on the dynamical mass. The most massive galaxies 
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have evolved passively at least since z ~ 1.2 while less massive systems show strong evidence 
of recent star formation activity, consistent with a “downsizing” scenario (Cowie et al. 1996). 
Within the range of environments covered by the GOODS-N held, no dependency is seen 
between this recent activity and the local galaxy density - contrary to the expectations of 
hierarchical models of galaxy formation (Diaferio et ah 2001; de Lucia et al. 2004). 

In this paper, we considerably extend and quantify the arguments presented in T05 in 
the context of a unifying hypothesis of the mass assembly history of spheroidals of various 
masses. We describe in detail our HST and Keck spectroscopic data and present the as¬ 
sociated catalogs. We derive the evolution of mass-to-light ratios based on analysis of the 
Fundamental Plane, taking into account selection effects and correlating deviations from the 
overall trend with morphological substructures. 

A plan follows. In §2 we introduce the GOODS-ACS multicolor imaging dataset, the 
construction of a photometric catalog and the procedure whereby which morphological types 
were assigned to all galaxies brighter than F850LP<22.5 (AB; hereafter zg). Surface photom¬ 
etry and scale-lengths are £t to the sample of spheroidals and morphological substructures 
identified via pixel-by-pixel color imaging (c.f. Menanteau et al. 2001). §3 presents the 
Keck spectroscopy and its reduction and the techniques used to measure stellar velocity 
dispersions. A complete catalog of the spectroscopic and photometric data is presented. §4 
describes the correlation of our catalog with the 2Ms X-ray catalog of the Chandra Deep 
Field North (Alexander et al. 2003; Barger et al. 2003) to identify nuclear activity. In §5 
we analyze the zero point, tilt, and scatter on the fundamental plane as a function of red- 
shift, and we use spectral and photometric diagnostics to provide independent evidence for 
recent star formation. §6 summarizes our findings and discusses our results in the context 
of formation scenarios for spheroidals. 

We assume throughout a cosmological model with Dm = 0.3, Da = 0.7 and Ho=65 
km s“^ Mpc“^. All magnitudes are in the AB system (Oke 1974), unless otherwise noted. 


2. Imaging 

Imaging data are taken from the public GOODS survey (Giavalisco et al. 2004). Ini¬ 
tial sample selection and planning of spectroscopic observations were performed before the 
completion of the survey and are based on the vO.5 release of epoch 1 (as described in 2.1). 
After completion of the GOODS-N survey and the release of the combined vl.O dataset, the 
photometric analysis was repeated using the newly processed and deeper images as described 
in this section. 
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2.1. Photometric catalog 

In preparation for our spectroscopic observations in April 2003, a first catalog was pro¬ 
duced by running SExtractor (Bertin & Arnouts 1996) on the vO.5 release of epoch 1 of the 
GOODS images. A zg selected catalog was constructed including all objects brighter than 
MAG_AUTO 2 : 9 ,vo .5 < 22.5. We will refer to this catalog as at the v0.5-zcatalog. After the re¬ 
lease of the combined GOODS images (release vl.O) a further catalog was constructed using 
the same procedure. We will refer to this catalog as the vl.O-zcatalog. The vl.O-zcatalog is 
available from Bundy, Ellis & Gonselice (2005; URL http: //www. astro. caltech. edu/G00DS_morphs/). 

Gomparisons between the two catalogs provide an estimate of the photometric accu¬ 
racy. After adjusting for the different zero points (c.f. the GOODS web page), the median 
difference and semi-interquartile range in zg magnitudes between vO.5 and vl.O are less than 
0.01 and 0.03 mags respectively. Since these are negligible with respect to other sources of 
errors throughout the paper, our catalog is effectively magnitude selected in vl.O as well. 
Accounting for the different zero points, the selection limit corresponds to Zg^^i.o < 22.43. 
Throughout the rest of the paper, all photometry is based on the vl.O data and zero points. 

Absolute magnitudes can be obtained from apparent magnitudes and colors using the trans¬ 
formations described in Appendix A. 


2.2. ACS Morphologies 

In view of the spectroscopic run scheduled in March 2003, an initial morphological clas- 
sihcation was undertaken by one of us (RSE) by visual inspection of the vO.5 epoch 1 2 : 9 -band 
images. The zg images were chosen to provide the closest match to optical rest frame at a 
redshift z ^ 1, minimizing the so-called morphological /c-correction. Morphologies (hereafter 
T) were assigned according to the Medium Deep Survey scheme introduced by Abraham 
et al. (1996a; see also Treu et al. 2003): T=-2=star, -l=compact, 0=E, 1=E/S0, 2=S0, 
3=Sa-|-b, 4=S, 5=Sc-|-d, 6=Irr, 7=Unclass, 8=Merger, 9=Fault. When the the deeper vl.O 
data release became available, a blind re-classihcation was also performed. Throughout the 
rest of the paper (unless otherwise specihed) and in the published catalog, all morpholog¬ 
ical types are based on the deeper vl.O data. Early-type and Sa-fb galaxies in the hnal 
spectroscopic sample are shown in Figures 2 and 3. 

Gomparison between the two morphological catalogs provides a check on the internal 
consistency of the classihcation and on its sensitivity to signal-to-noise (Figure 1). Gonsid- 
ering the entire sample, the agreement is very good, with a median offset of 0, a mean offset 
of -0.1 and a standard deviation of 1.3 (in the 12 point T scheme described above). However, 
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for the class of spheroidal objects with T=0,l,2 in the vO.5 catalog, the mean offset increases 
to 0.4 (median=0, standard deviation 1.2). The differences arise because of the improve¬ 
ment in signal-to-noise in the vl.O images which reveal, when present, the characteristic low 
surface brightness features of spiral disks. Thus objects classed as spheroidals in the vO.5 
catalog (which formed the basis of our spectroscopic sample) include some systems which 
later turned out to be bulge-dominated spirals, as described quantitatively in Section 3.1. 
Of course, by including these outliers in our sample, we have a valuable measure of the 
uncertainties associated with dehning our morphological boundaries. 

Although we include our morphological classihcations in the large catalog discussed by 
Bundy, Ellis & Conselice (2005), it is important to examine the extent to which classihcation 
uncertainties permitted in large statistical analyses might influence the results of this, more 
focused, study of a relatively small number of spheroidals. A particular concern is the 
likelihood of contamination by the more numerous bulge-dominated spirals. 

We can estimate the probability of spiral contamination as follows. At the magnitude 
limit of the survey, the rms scatter rmstot in classihcation between the vO.5 and vl.O catalogs 
is 1.3 types. Assuming that uncertainties scale as the square root of the exposure time, we 
obtain rmstot=\/rmsf q -|- rmsg 5 =A/rmsf q -|- 5rmsf q, since the vl.O images represent 5 times 
the integration of the vO.5 images. This would appear to be a conservative assumption, 
since the vl.O images also have improved sampling, defect removal and relative alignment. 
Thus, we expect rmsi.o=rmso. 5 /y^O) ~ 0.5 implying that ~l/6 of the Sa-l-b sample might 
be misclassihed as SOs. As Sa-|-bs are twice as numerous as SOs overall, we estimate that 
~l/3 of the SOs could be misclassihed Sa-|-bs, i.e. approximately 20. 

To identify possible misclassihed Sa-|-bs, residuals from the htting (faint analogs 
of the ones shown in Figure 4) were visually inspected by two of us (RSE,TT) for disk 
asymmetries and spiral features. To be conservative, all edge-on SOs, where inclination could 
be hiding spiral features, were also hagged. With these criteria, we identihed 23 possible 
Sa+b contaminants (ID 199, 312, 369, 400, 402, 424, 505, 63, 655, 786, 811, 870, 872, 895, 
996, 1286, 1441, 1463, 1485, 1491, 1577, 1685, 1764). The technique is surprisingly ehective 
in locating faint spiral features but we use this bagging only to illustrate the ehect of possible 
contamination rather than to rehne the published morphological catalog. We also identihed 
7 galaxies with prominent dust lanes or asymmetries (ID 312, 377, 460, 879, 1354, 1477, 
1764) where the dust and asymmetries could bias photometric and spectroscopic parameters 
used in our analysis. In Section 5.3 we will test the robustness of our analysis with respect by 
these possible biases by comparing results obtained with and without these two subsamples. 
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2.3. Surface photometry 

Surface photometry of all early-type galaxies and bulge dominated spirals targeted for 
spectroscopy was performed using the publicly available software GALFIT (Peng et ah 2002). 
Effective radius (Re) and magnitudes (and hence effective surface brightness SBg) were ob¬ 
tained in all four bands as the best fitting parameters of a surface brightness profile, 
following the standard procedure used to construct the Fundmental Plane (Jprgensen, Franx 
& Kjaergaard 1992; van Dokkum & Franx 1996; Treu et ah 2001). The fitting procedure 
requires a point spread function (PSF) to convolve the model before computing the . For 
each galaxy we adopted the nearest bright non-saturated star in the GOODS-N field as 
fiducial PSF, thus correctly taking into account the modifications of the PSF due to the 
GOODS vl.O reduction procedure. Examples of fits and residuals from a profile are 

shown in Figure 4. 

As is commonly known (e.g. Kelson et al. 2000a; Treu et al. 2001a) the particular 
combination of photometric parameters necessary for constructing the FP: FPph = logi?e ~ 
0.32S'Re, is limited for datasets of good quality such as GOODS by the fitting technique. 
Tests with fit2D (Treu et al. 2001a) indicate that the FPph parameter derived using GALFIT 
is consistent with those derived in previous studies to a level of 0.02 or better. An additional 
test can be obtained by comparing our photometry to the one of vDE03, for the 9 objects 
in common. The effective radii agree extremely well (average difference log Re=0.0004, rms 
scatter 0.066), and the photometry is in good agreement (F814W-F775=0.010 mag, with 
rms=0.34 mag) considering the different filters and instruments (assuming F814W=F775W). 
As a result FPph agrees well (average difference 0.019, rms 0.051). 

We adopt 0.03 as our best estimate of the systematic uncertainty on the FPph due to 
the fitting technique. Absolute magnitudes can be obtained from apparent magnitudes and 
colors using the transformations described in Appendix A (c.f. Treu et al. 2001a). Rest 
frame effective radii are obtained by interpolation from the two nearest filters as described 
in Treu et al. 2001a. 


2.4. Internal colors, bine cores and morphological asymmetry 

The GOODS images are, like their HDF predecessors, of adequate signal to noise for us 
to examine the internal homogeneity as a verification of the morphological classifications and 
as a possible diagnostic of recent activity. Of particular interest, given our large sample with 
associated spectroscopy are the so-called ‘blue-core’ spheroidals first identified by Menanteau 
et al. (2001) in the HDF images. In this paper we define a more robust way to uniformly 



select such systems at various redshifts and our resulting sample is significantly larger than 
any hitherto. 

Rest frame pixel-by-pixel color maps (B-V) were derived by applying the transformations 
described in Appendix A, with the aim of identifying internal color peculiarities which might 
indicate recent starformation activity. The vast majority of early-type galaxies proved to 
be relatively smooth and homogenous. However 14/163 early-type galaxies with measured 
redshift show a clear blue core. 

Menanteau et al. (2001) already addressed the key issue of whether the blue cores could 
be artefacts arising from the wavelength dependency of the HST PSF. Although that analysis 
was conducted on the basis of WFPC-2 images in the Hubble Deep Field and rich clusters, 
the improved spatial sampling of the ACS makes their arguments tighter in this case. An 
instrumental origin would imply a much greater frequency of blue cores at a given apparent 
magnitude and not produce the environmental dependencies discussed by Menanteau et al. 
or the correlations with independent diagnostics of star formation discussed below. 

We chose to define a blue core spheroidal as one where a difference of at least (5(B- 
V)=0.2 is seen between the core of the galaxy and the outer parts. As our color selection is 
made in rest-frame color, this is significantly more useful than the observed V — I criterion 
used in previous studies (Menanteau et al. 2001, 2004, 2005). We find a total of 14 blue-core 
early-type galaxies in GOODS-N and these are sorted by redshift and shown in Figure 5. 
The nature of these systems is discussed further in Sections 4 and 5.6. 


3. Keck Spectroscopy 
3.1. Sample selection and mask design 

The sample for spectroscopic follow-up was selected from the v0.5-zcatalog according to 
the following criteria: ; 29 ,«o .5 < 22.5, —1 < T < 2. Self-consistent astrometry across the field 
was obtained by matching the ACS catalog to the ground based catalog used by the Team 
Keck Redshift Survey (TKRS; Wirth et al. 2004). Three sets of masks with their major axis 
aligned with that of the GOODS-N mosaic were designed (north, center, south). For each 
pointing a first submask was filled with objects at random from the sample of early-type 
galaxies using the vO.5 morphologies. A second submask was obtained by replacing objects 
brighter than 2 : 9 ,„o .5 = 21.5 with new objects. The goal of this procedure was to obtain longer 
exposure times for the faintest objects while maximizing the number of survey objects. 

In total, 261 galaxies from the primary spectroscopic sample were included in the six 
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masks, plus 22 T=3 (vO.5) galaxies. Based on the subsequently-determined deeper vl.O 
morphologies, the spectroscopic sample 175 early-types, 66 early-spirals (T=3), 14 later 
type galaxies (T>4) and 28 stars. This mixing of types is expected because of the increased 
signal-to-noise in the vl.O (see section 2.2) images and because the number of objects that 
can drift out of the sample due to random errors in the classihcation is much larger than 
the number of objects that can enter, and because a conservative star/galaxy separation 
criterion was used in the vO.5 images. 

Eighteen slitlets (with a random distribution of morphologies and luminosities) were not 
milled or could not be extracted because they fell too close to the edges of the masks, so the 
hnal tally is comprised of 265 objects. Six objects were observed twice, in different slitlets 
of different masks, to allow for internal checks. All the galaxies with measured redshifts are 
shown in Figure 2 and Figure 3. 

To ensure that the retrospective vl.O morphological cut does not bias the properties 
of our sample we compare in Figure 6 the luminosity and color distribution of the parent 
sample comprising all early-type galaxies in the GOODS-N field (vl.O, empty histogram) with 
dehned in onr spectroscopic sample (vO.5, hatched histogram). The color distribntions are 
very similar with a typical sampling rate between 50 and 70%, showing that onr spectroscopic 
sample is nnbiased in color. However, the sampling rate is a declining fnnction of apparent 
magnitnde beyond zg ~ 21.5, as expected, because of our mask-design strategy, the higher 
snrface density and the higher probability of morphological misclassihcation. 

In the following we will consider the early-type galaxies as onr primary sample. The 
early-spiral sample (composed mostly of objects classihed as early-type in the vO.5 analysis) 
will be useful as an ancillary in order to compare with results based on shallower imaging 
data. A complete discussion of the properties and star formation histories of the bnlges will 
be given elsewhere. 


3.2. Observations and data reduction 

The held was observed using the Deep Imaging and Mnlti-Object Spectrograph (DEIMOS) 
at the Keck-II telescope on the nights of April 1 to 5 2003. Conditions were excellent, with 
clear sky and seeing in the range Cf.'O-O'.'S thronghont the rnn. The 1200 grating blazed at 
7500A was centered at 8000 A providing a pixel scale of 071185 x 0.33A, and a typical 
wavelength coverage of 2600 A. The 170 wide slitlets yielded an instrumental resolntion of 
(Ts ~ 20 — 30 kms“^, as determined from arc lines and sky emission lines. Total exposure 
times ranged between 14,400 and 38,900 seconds as listed in Table 1. 
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The data were reduced using the IDL pipeline developed by the DEEP2 collaboration 
(Davis et ah 2003; Faber et ah 2005, in prep). After extraction, the Id spectra were assigned 
a redshift using the software developed by the DEEP2 collaboration (described in Coil et ah 
2004). All spectra were examined independently by two of us (TT and TXL) and a quality 
criterion was assigned (q=-l stars, q=2 possible, q=3 likely, q=4 certain). In those very few 
cases where the redshift/quality parameter differed, a simultaneous inspection by the two 
classihers was sufficient to reach an agreement. 

Redshifts were obtained for all the 265 targeted objects, including 217 from the primary 
sample, 21 Sa+b galaxies, and 27 that turned out to be stars (vO.5 classihcation). Based 
on the superior vl.O classihcations, the sample with redshifts is comprised of 163 early- 
type galaxies, 61 early-spirals, 23 later type galaxies and 1 star (according reclassihed as 
compact). Comparison with other surveys shows that our redshifts are accurate to 30 kms“^ 
or better (c.f. Wirth et al. 2004). An approximate relative flux calibration was obtained 
using the measured sensitivity curves for DEIMOS. Examples of spectra covering the range 
of redshifts and luminosities (and hence signal to noise ratio) are given in Figure 7. The 
redshift distribution of the E-I-SO and Sa-l-b samples is shown in Figure 8. 


3.3. Internal kinematics 

Stellar velocity dispersion were obtained by comparing the extracted spectra with spec¬ 
tra of G-K giants and of the Sun (a high resolution solar spectrum was obtained from the 
URL http://bass2000.obspm.fr/solar_spect.php) as described in detail in Treu et al. 
(1999; 2001) and in Sand et al. (2004). A brief summary follows. 

First, the instrumental resolution was measured by htting gaussians to well-exposed and 
unblended sky lines and arc calibration lamps. This corresponds to a velocity dispersion of 
(Js = 30 — 20 km s“^in the spectral range of interest, and is well described by a second order 
polynomial in wavelength (Js(A). Second, high resolution template stars were redshifted and 
smoothed to match the resolution of the instrumental setup, and convolved with gaussians 
in log A space to reproduce the kinematic broadening. The best fitting velocity dispersion 
was found using the pixel htting code developed by van der Marel (1994). Spectral regions 
affected by emission lines, Balmer lines, atmospheric A and B band absorption, as well 
as residual defects were carefullly masked out during the ht. Those pixels affected by the 
strongest sky emission lines and their immediate neighbors were similarly hagged. Each 
object was htted independently with each template and using different spectral regions. The 
G-K giant library extends from 4000 A redwards, while the solar spectrum extends to the 
UV atmospheric limit, enabling us to take advantage of near-UV absorption features in the 
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highest redshift objects (c.f. van Dokkum & Stanford 2003). 

Each spectrnm and htting result was carefully inspected by eye by one of us (TT) 
to identify and correct potential problems (such as insufficient wavelength range or reduc¬ 
tion/extraction defects). Two example of spectra with kinematic fits - an intermediate 
quality spectrum and the highest redshift object with measured velocity dispersion - are 
shown in Figure 9. 

The majority (181/224 of all galaxies; 141/163 spheroidals) of objects yielded a stable 
velocity dispersion measure ((Tap) listed in Table 1. A few objects failed to deliver a reliable 
velocity dispersion, due to insufficient signal to noise or unfortunate location of the strongest 
absorption features on persistent sky subtraction residuals. For a handful of objects at 
low redshift the procedure yielded velocity dispersions below the instrumental resolution 
(Ts = 30 km s“^. These measures were considered insufficiently reliable enough and discarded 
(although plausibly correct given their faint intrinsic luminosity and the fact that the NaD 
line is resolved). The success rate as a function of redshift, luminosity, and average signal to 
noise ratio (per A, observer frame) is summarized in Figure 10. For simplicity, we report the 
average signal to noise ratio across the entire spectrum as a measure of the total information, 
which is approximately given by <S/N> multiplied by the square root of the number of A, 
i.e. V2600. The <S/N> per A in the rest frame is readily obtained as <S/N>rest=<S/N> 
Vl + As other authors (e.g. Jprgensen, Franx & Kjaergaard 1995) often quote values 
based on smaller wavelength regions, or ones longward of the 4000 A break, in comparing 
the limiting S/N for which dispersions were successfully measured, differences in dehnition 
should be borne in mind (see Figure 11). 

For the FP analysis, the central velocity dispersions (i.e. within a circular aperture of 
radius rg/S) was obtained as (T=1.10±0.04 (Tap (Treu et al. 2001). 


3.4. Tests of the Measured Velocity Dispersions 

Given the size of our sample, a key issue before the scientific analysis that follows is the 
reliability of our stellar velocity dispersions. We have approached this via several independent 
tests which, collectively, demonstrate a high degree of precision has been attained, thanks 
largely to the long integrations and the excellent quality of the data obtained with the 
DEIMOS spectrograph. 

1. Repeated measures of six objects observed with DEIMOS in different slitlets/conhgurations 
(Figure 11) show that our measurements are stable with respect to internal effects such 
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as varying the wavelength range and the slitlet position in the held of view of the spec¬ 
trograph. 

2. A comparison with the 7 objects in common with van Dokkum & Ellis (2003) examined 
independently with the LRIS spectrograph shows that both sets are also statistically 
indistinguishable, with < daja >= —0.06 and an rms scatter of 12 %, consistent with 
the estimated errors (0.12/y/(6) = 0.05). 

3. A subsample of spectra were analyzed independently by one of us (PGD) using methods 
described in van Dokkum & Franx (1996). The internal comparison indicates stellar 
velocity dispersions that agree to better than 1% on average, with an rms scatter of 
12%, indicating that there are no signihcant systematic uncertainties due to the choice 
of htting method. 

We conclude that our estimated errors on velocity dispersion (on average less than 10%) are 
accurate. 


3.5. Basic Star Formation Diagnostics: [O II] and H(5 

Whenever the redshift placed the useful diagnostic features of [O II] 3727 A and H(5 
4101 A in the observed spectral region, we measured their equivalent widths using an IDL 
procedure. The band dehnitions of Fisher et ah (1998; [O II]) and Trager et al. (1998; H^a) 
were used. The noise properties of the spectra were taken into account so as to maximize 
the signal-to-noise ratio. Observed frame equivalent widths in A are listed in Table 1 In the 
following analysis we will always refer to rest frame equivalent widths in A. 


4. X-ray properties 

To obtain a complementary source of information on nuclear activity and star formation 
in our sample of spheroidals (and in particular on the blue cores) we cross correlated our 
redshift catalog with the published X-ray 503 sources catalog of the Chandra Deep Field 
North (Alexander et al. 2003; Barger et al. 2003), Ending 42 sources with X-ray counterparts 
(26 E-|-S0s; 13 Sa-|-bs; 3 Ss). All the 42 sources have redshifts published in Barger et al. 
(2003), in agreement with our determinations. 

From the soft and hard fluxes measured in Barger et al. (2003), we derive the rest 
frame soft (0.5-2 keV) and hard (2-8 keV) emitted luminosities applying a /c-correction for 
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a spectrum with photon index F = 2. A spectral slope of F = 2 is appropriate for X- 
ray emission from prompt star formation. Since the majority of the 42 X-ray counterparts 
of our spheroidal galaxies have emitted luminosities in the range typical of star forming 
galaxies (Fo. 5_2 < 10"^^ erg s“^, see Ranalli et ah 2003), we concluded that F = 2 is a 
good approximation of the required /c-correction. Adopting a different /c-correction has a 
negligible impact on our results, as described in the next paragraph. 

For simplicity, we will classify our objects as AGN-dominated if the X-ray luminosity 
(soft or hard) is above 10^^ erg s“^, and as star formation dominated or low luminosity 
AGN (LLAGN) if it is not in the previous category and the X-ray luminosity is above 
10^0 gj,g g-i following we will label these two groups as XAGN and XSF/LLAGN; 

motivations for these working dehnitions is given in Appendix B). With these dehnitions, 17 
sources (9 E-|-S0s) are XAGNs while 21 (16 E-|-S0s) are XSF/LLAGNs. Adopting different 
slopes for the /c-correction, e.g. F ~ 1.4 typical of the average spectrum of absorbed and 
unabsorbed AGN that dominate the population in the faint flux regime (as in the GFDN), 
would only adjust these numbers in a minor way to 16 and 20 galaxies respectively in the 
XAGN and XSF-LLAGN categories. The same classihcation is obtained by adopting the 
spectral slope evaluated for each source from the X-ray hardness ratio (see Alexander et al. 
2003). Interestingly, only 2 of the 14 blue cores E-|-S0s are classihed as XAGN (ID=1323,872 
Figure 5), while 1 is classihed as XSF-LLAGN (ID=1362), showing that most blue cores are 
not due to a central AGN. 

For the X-ray sources classihed as XSF/LLAGN we can obtain an upper limit to the star 
formation rate using the conversion M = 2.2 x 10“^‘^Lo.5-2fcey Mq yr“^ (Ranalli, Gomastri 
& Setti 2003). The upper limits are in the range 1-100 Mq yr“^ with an average (median) 
value of 38±30 Mq yr“^ (28±23 Mq yr“^). We would like to emphasize that in many of 
these cases, these upper limits will be well above the actual star formation rate, because of 
the LLAGN contribution and/or emission from X-ray halos and LMXB. This is conhrmed 
also by the optical spectra which do not show the dramatic features that would be expected 
for star formation rates of order ~20 Mq yr“^. In order to disentangle more accurately the 
various contributions - especially the soft extended halos from the typically hard and point- 
like AGN - a detailed analysis of the spectrum and the spatial distribution of the X-ray 
emission will be necessary. This goes beyond the scope of the present paper and is left for 
future work. 
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5. Results 

In this section we use our measurements to study the evolution of the stellar populations 
of early-type galaxies, using the FP as a primary diagnostic tool. In § 5.1 we derive the 
formalism to intrepret the evolution of slopes, intercept and scatter of the FP in terms of 
(mass dependent) evolution of the effective mass-to-light ratio. In the following sections, 
§ 5.2 and § 5.3, we will use the FP to study, hrst qualitatively and then quantitatively, 
the evolution of the mass-to-light ratio of spheroidal galaxies, showing that both the slope 
and scatter of the FP change with redshift. We will interpret these results in terms of the 
evolution of the stellar populations. In § 5.5 we will compare our results with the existing 
literature, showing that previously considered discrepant results can be reconciled if the 
proper selection effects and uncertainties are taken into account. Finally, in § 5.6, we will 
examine the relationship between the evolution of the mass-to-light ratio and other stellar 
population diagnostics, providing independent evidence for recent starformation. 


5.1. The Fundamental Plane as a diagnostic of stellar populations 

The Fundamental Plane is dehned as 


logi?e = a log a + (3 SBe + 7 , 


( 1 ) 


(a in km s“\ Re in kpc, SBe in mag arcsec”^). For this study we will adopt as the local 
relationship the FP of the Coma cluster {a = 1.25, (3 = 0.32, 7 = —8.970, in B(AB) for 
Ho = 65km s“^Mpc“^; Jprgensen, Franx & Kjaergaard 1996, hereafter JFK96). We use Coma 
as the local reference for both cluster and held to minimize systematic uncertainties related 
to hlter transformations, distance determination, and selection effects; see discussion in Treu 
et al. (2001b). However, the effect of the systematic uncertainty in the local intercept will 
be discussed and quantihed in the analysis. We will also adopt 0.08 in logi?e as the intrinsic 
scatter of the local relationship, consistent with cluster and held estimates (JFK96, Bernardi 
et al. 2003). 


A simple physical interpretation of the FP (Faber et al. 1987) can be given by dehning 
an ehective mass, 

, ^ 5(7^i?e / \ 

M = ( 2 ) 

(c.f. Bender, Burstein & Faber 1992), and by dehning the luminosity in the usual way 


—2.5 log L = SBe — 5 log Rq — 2.5 log 27r, 


( 3 ) 
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and by considering an effective mass-lnminosity relation of the form (also known as the “tilt” 
of the FP; e.g., Ciotti, Lanzoni & Renzini 1996; Lanzoni et ah 2004) 

L oc M^. (4) 

These assnmptions lead directly to the FP relation given above, provided 

a-10/3 +2 = 0, (5) 

(van Albada, Berlin & Stiavelli 1995), with rj = 0.2q;// 3, which is very closed to the observed 
valne. When converted in solar nnits and in the form of eqnation 4, reads: 

log M/Lb = 0.25 log M - 1.934. ( 6 ) 


In this framework, variations of the slopes a, (3 and the intercept 7 as a fnnction 
of redshift are easily interpreted as general trends in Inminosity evolntion of the stellar 
popnlations. If a and Re do not evolve with redshift, for an individnal galaxy (labeled by 
the snperscript i) 

y = logR(,-aloga*-/3SB(,, (7) 


the offset with respect to the prediction of the FP (Ay* = 
the M/L by 


A log 



Ay* 


y* — y) is related to the offset of 


( 8 ) 


which can be used to measure the average evolution and/or scatter of M/L at given M (Note 
that equations 7 and 8 are independent of the evolution of slopes and that this relation does 
not require the FP to be of the form in Equation 6 .). Therefore, the tightness of the FP 
(JFK96; Bernardi et ah 2003) constrains the homogeneity of the stellar populations of early- 
type galaxies (similarly to the small scatter of the color-magnitude relation. Bower, Lucey 
& Ellis 1992). An evolution with redshifts of the slopes of FP can be interpreted in terms of 
a mass dependent evolution of M/L. For a FP of the form in Equation 6 the evolutionary 
rate depends only on effective mass: 


d\og{M/L) _ 1 

In the following, we will often compare the evolution of the FP for field early-type 
galaxies to that found for cluster E+SOs. For this purpose, we will adopt as our fiducial 
evolution for E+SO in massive clusters 


da y — log M — log G + log 5 dy 
dz 10/3 dz 


( 9 ) 


d\og{M/L- q)/ dz = —0.46 ± 0.04 


( 10 ) 
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(van Dokkum & Stanford 2003; see also van Dokknm et al. 1998, who find —0.49±0.05). We 
also define a parameter that measnres the differential evolntion between the mass to light 
ratio of a held E±S0 and that expected for the hdncial cluster E±S0: 

5A log(M/LB) = A log(M/LB) ± 0.46^. (11) 


5.2. Evolution of the FP: qualitative results 

The location of the E/SO (pentagons) and Sa±b (stars) galaxies in the FP-space is 
shown in the hve panels of Figure 13, each panel showing galaxies in a different redshift bin. 
The local relationship measured in the Coma cluster (JFK96) is shown as a solid line for 
reference (see Section 5.1 for discussion). The early-type galaxies dehne a relatively narrow 
FP in each redshift bin (note that the hnite size of the bin introduces some artihcal scatter), 
while the scatter is larger for Sa±b galaxies. As illustrated in T05, galaxies move away from 
the local relationship as redshift increases, becoming brighter at hxed effective radius and 
velocity dispersion. 

The offset from the local relationship can be converted into a change in effective mass-to- 
light ratio using Eqn 7 and 8. The change for each individual galaxy as a function of redshift 
is shown in Figure 14. Evolutionary tracks for a simple (single-burst) stellar population ^ 
formed at Zf = 1, 2, 5 are also shown for comparison. It is clear that even the genuine E±S0 
population spans a wide range in Ay^, covering the entire range predicted by these simple 
models. We show these models only for illustration at this point; such formation epochs 
would refer approximately to the luminosity weighted age of the stellar population. Even a 
modest fraction of younger stars would have a signihcant influence on the mass-to-light ratio 
(Trager et al. 2000; Treu et al. 2001b). 

As discussed in T05, the observed change in the slope in Figure 13 and increased scatter 
in Figure 14 with respect to local samples, are the result of mass-dependent trends, combined 
with selection effects and increased intrinsic scatter of the FP. This is illustrated in Figure 14 
by plotting the offset from the local FP as a function of galaxy mass. The most massive 
galaxies (M> IO^^'^Mq) have relatively high mass-to-light ratios, minimal scatter and follow 
the passively-evolving track with a high formation epoch fairly closely . As less massive 
galaxies are included, the scatter and the average offset from the local relation increase. 
This possible trend was already tentatively seen in previous studies (Treu et al. 2002; van 


^Unless otherwise stated all evolutionary models are computed using Bruzual & Chariot 1996 stellar 
population synthesis models GISSEL96, Salpeter IMF, solar metallicity and Kurucz atmosphere models. 
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der Wei et al. 2004; see Pahre 1998, Kelson et al. 2000c, Wuyts et al. 2004, Moran et al. 
2004 for similar results in the cluster context), but only with our large sample does it seem 
conclusive (T05). In the next section we will obtain quantitative limits on the evolution of 
the mass-to-light ratio and its scatter as function of mass and redshift taking due account 
of selection biases. 


5.3. Evolution of the FP: quantitative results 

As we have seen in the previous section, the evolution of the FP depends on galaxy 
mass, therefore we need to allow the FP slopes to vary with redshift. In addition, since at 
face value it appears that the intrinsic scatter around the FP increases with redshift, the 
question of possible evolution in the intrinsic scatter is of great interest. As shown in Treu et 
al. (2001b; 2002) the evolution of intercept, slope and scatter are interconnected and deeply 
dependent on selection effects. The larger the intrinsic scatter, the larger the Malmquist- 
like bias on the evolution of the mass-to-light ratio, and therefore the larger the correction 
necessary to recover the “unbiased” evolutionary trend. 

The bias due to luminosity selection is illustrated in Figure 15 where we show the 
location of the early-type galaxies in a different projection of the FP: the M — M/L plane 
(analogous to the k space introduced by Bender, Burstein & Faber 1992). The dotted line 
represents the local relationship, while as in Figure 13 the pentagons represent the early-type 
galaxies. Hatched regions represent those excluded by our magnitude limit. As expected, 
objects with low mass and high mass to light-ratio are excluded a priori from the sample. 
Only the two lowest redshift bins probe deep enough to sample the entire FP. Ai z ^ 0.5 and 
above, selection effects begin to dominate the low mass trends. Although we now embark 
on a rigorous modeling of the biases, we can already infer from the existence of galaxies 
brighter than the selection limit that the intrinsic scatter probably increases with redshift 
and therefore that the smaller mass galaxies have a diverse star formation history. 

We note, in passing, that it is not sufficient to deduce evolutionary trends by considering 
the high-mass objects only (e.g. van der Wei et al. 2004), since the uncertainties on M/L and 
M are correlated. Selecting objects above a certain mass would introduce an additional bias, 
leading in general to an overestimate of the M/L and hence an underestimate of the evolution 
of this subset of the population. In short, there is no way to avoid a proper modeling of the 
observables. 



5.3.1. A Bayes an-Montecarlo approach 


Selection effects can be properly taken into account by using the Bayesan-Montecarlo 
method introduced by Treu (2001) and Treu et ah (2001, 2002). Briefly, the method uses 
a Montecarlo algorithm to compute the likelihood of obtaining the observed surface bright¬ 
nesses for hxed velocity dispersion and effective radius, given the uncertainties and sample 
selection properties as a function of the FP coefficients and scatter. The likelihood is then 
combined with a prior using Bayes Theorem (1768) to recover the posterior probability for 
each parameter. The likelihood is generally non-zero over a small region of the possible 
parameter space and therefore the precise choice of the prior is not critical, assuming that 
is chosen to be regular enough. In the following, we adopt a uniform prior (see Treu et al. 
2001 for a discussion). All results presented in this refer only to the spectroscopic sample of 
E-I-SO galaxies, that is spirals are excluded. 

In this section we will derive formal limits on the evolution of the coefficients of the FP 
using increasingly complicated models. As a hrst step, we consider a simplihed model where 
only the FP intercept 7 and scatter are allowed to vary linearly with redshift, and we thus 
seek limits on the derivative. This is the traditional approach taken by most studies and 
therefore it is useful for comparison with previous work. This procedure gives an estimate of 
the average mass to light ratio evolution of the sample and an upper limit to the evolution 
of the scatter. By construction the scatter for hxed slopes is always larger than the scatter 
around the best htting slopes. The posterior probability contours for our sample in this case 
are shown in Figure 16. 

As expected the constraints on the two quantities are correlated in the sense that less 
evolution in the mass to light ratio is required for a population with larger scatter. The 
hypothesis that the FP evolves with a hxed slope and constant scatter is rejected by our 
dataset. The best estimate and 68 % limits (obtained as described in TOl) are dj/dz = 
0.58^o;og and da.y/dz = 0 . 18 ^q;o 4 corresponding to < d\og{M/L-^)/dz >= —0.721q;q 5 and 
duiogM/i/dz < 0.26. As shown in Figure 16, the average evolution of the mass-to-light ratio 
for our sample of held E-|-S0s is signihcantly faster than that of cluster E-|-S0s. 

To verify that our results are not inhuenced by the inclusion of galaxies with dust- 
lanes or possible early-type spirals, we re-ran our Montecarlo simulations excluding the 
objects hagged in Section 2.2. Excluding the dusty/peculiar objects (5 with measured 
velocity dispersion), changes the posterior probability in a negligible way. Excluding the 
possible Sa-|-b galaxies (22 with measured velocity dispersion) changes the above limits to 
< d\og{M/LB)/dz >= -0.75^°;^° and daiogM/i/dz < 0.28, i.e. much less than the uncer¬ 
tainties (see Figure 16). We conclude that uncertainties in the morphological classihcation, 
although likely at the level discussed, do not ahect our results. 



In addition, we compute the effects of systematic uncertainties on the local intercept. 
For simplicity, we will express variations in the local intercept as variations in surface bright¬ 
ness SBg at a hxed velocity dispersion and effective radius (c.f. Bernardi et al. 2005). Since 
the evolution depends linearly on the local intercept, there is a simple linear relation be¬ 
tween uncertainties, i.e. 5dlog{M/ L-q) /dz = -0.57 (5SBe, where the coefficient is a weighted 
average of individual redshifts and offsets (the coefficient has been verihed numerically by 
repeating the Montecarlo simulations after varying the local intercept). The environmental 
dependence of the intercept in the local universe is very small and hard to quantify because 
of selection effects, hlters and parameters dehnitions. As a conservative approach, we will 
adopt the maximum environmental difference observed in the Sloan data (Bernardi et al. 
2005), corresponding to (5SBe = 0.075 mags, i.e. a systematic error on the evolutionary rate 
of ±0.04. 

Next, we consider a more realistic evolutionary scenario: namely evolution of the slopes 
a and (3, and hence mass dependent evolution (Equation 9). In implementing this model, 
we assume that the FP can be described with a simple power law relation between M and 
M/L at any redshift, i.e. a{z) — 10/3(2:) ±2 = 0. We will also assume that evolution of the 
coefficients and slopes can be described by a linear relation, e.g. a{z) = q;(0) ± da/dz\QZ. 
Thus our model has three free parameters, da/dz\o, dj/dzlo and da^ / dz\o. For simplicity, 
we will drop the index 0, assuming all derivatives are computed at z = 0. 

Posterior probability contours (68% and 95%) for da/dz and d'y/dz, marginalized over 
da^/dz, are shown in Figure 17. Constant slopes are clearly rejected, demonstrating unam¬ 
biguously that the tilt of the FP evolves with redshift. This result is translated in terms 
of effective mass to light and mass in Figure 18 using Eqn 9. The evolution of the mass to 
light ratio depends strongly on mass: the most massive galaxies evolving similarly to massive 
cluster galaxies (i.e. d\og{M /L-q)/ dz = —0.46), while evolution is faster for smaller masses. 
This is the downsizing trend discussed in T05. 

As a hnal step, we derive the posterior probability distribution function for the evolution 
of the intrinsic scatter of the FP da^/dz - marginalized over da/dz and d'y/dz - shown in 
Figure 19. The peak value of the probability distribution of function and the 68% contours 
are da^/dz = 0.032 ± 0.012, corresponding to dimsilog M/L-q)/ dz = 0.040 ± 0.015. Note 
this measurement of the scatter is much smaller than the scatter at fixed slopes (Figure 16), 
showing that part of the visual impression of the scatter when the evolution is plotted as in 
Figure 14 is dominated by the change in the slope. This proves that the scatter in the FP, 
and therefore the scatter in stellar population properties, increases with redshift, as expected 
for example if the scatter is due to a spread in stellar ages. 
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5.4. Implications for the star formation history of E+SO galaxies 

In this section we proceed to interpret our findings on the evolution of the FP in terms 
of star formation and assembly history of E+SO galaxies. Since only objects morphologically 
identified as E+SO are included in our sample, what we describe is the effective star formation 
history of the existing population at any given redshift. If there was significant morphological 
evolution (“progenitor bias”, as discussed by van Dokkum & Franx 2001 for cluster early- 
type galaxies), the average star formation activity of the parent population would be more 
delayed. 

First, we will use single burst stellar population models showing that they imply very 
low redshift of formations for the population with effective mass below 10 ^^ Mq, inconsistent 
with independent star formation diagnostics and other observational facts. Then, we will 
proceed to explore more realistic models of composite stellar populations. 

The traditional physical description of the evolution seen in the FP is the change in mass 
to light ratio expected in the redshift range z = 0 — 1.2 for a simple (single-burst) stellar 
population (SSP) (the correspondence between d\og{M/L b)/ dz and redshift of formation 
Zf is shown in Figure 20). For early-type galaxies more massive than the evolution 

implies a high redshift of formation, Zf > 2, whereas redshifts of formation as low a.s Zf = 
1.2 — 1.3 are adequate for the systems of lower effective mass. 

A major result of our paper, however, is that the above description is too simplistic 
because various independent diagnostics available to us (see Section 5.6) indicate that a 
significant fraction of the early-type population with mass lower than 10 ^^ Mq has suffered 
recent/secondary star formation. Furthermore, a SSP redshift of formation of Zf ~ 1.2 for 
lower mass objects would be impossible to reconcile with the results of deep infrared surveys 
(Glazebrook et al. 2004; Cimatti et ah 2004), which have found a substantial number of 
galaxies (including morphological E+SO; e.g., Stanford et al. 2004) with stellar mass above 
10 ^° Mq at redshift above 1.2. 

A more appropriate representation is therefore a composite stellar population. Since 
subsequent bursts involving a relatively small fraction of stellar mass can substantially alter 
the evolution of the mass to light ratio (Trager et al. 2000; Treu et al. 2001), we will consider 
models where a (dominant) fraction of the stellar mass is old (zfi > 2) and the rest of 
the stellar mass was formed at Zf 2 < 1.2. We will show that these models provide a good 
description of the data and use them to derive limits on the fraction of younger stars. 

This procedure is illustrated in Figure 21 where we compare the evolution of the mass- 
to-light ratio for families of simple two-burst models with that observed for our sample. 
Starting from the left panel panel we discuss E+SO galaxies in order of decreasing mass. For 
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the most massive galaxies (> the average evolution is well described by that of a 

very old stellar population formed at Zf = 5. Secondary burst at Zf 2 < 1.2 are allowed but 
they cannot involve more than 1% of the total stellar mass. For intermediate mass objects 
(center panel), larger bursts are allowed, up to ~ 5 % of the stellar mass. For the smallest 
masses (right panel) bursts involving 10% of the stellar mass at z < 1.2 are needed. For the 
largest masses these are solid upper limits and are consistent with, for example, accretion 
of small satellites. As discussed in T05 this is at variance with the expectations of a simple 
scheme where the growth of stellar mass follows the bottom-up growth of dark matter halos 
(see T05 for discussion). 

For the less massive galaxies, however, when the mass in the secondary bursts starts to 
be a signihcant fraction of the total mass there is an additional complication. These larger 
secondary bursts are a quite dramatic event, which could be the result for example of a major 
merger, involving a rapid transformation of gas into stars. After such a dramatic event, it 
is unlikely that a galaxy would be recognizable morphologically until dynamical relaxation 
is complete (i.e. of order a Gyr after the event; c.f. van Dokkum & Franx 2001). Thus a 
more realistic estimate of the mass in these signihcant secondary bursts is obtained as in 
T05, assuming at least IGyr of age for the young population. In some cases, this requires 
fraction of mass in the young population as high as 20-40%. 


5.5. Comparison with previous results 

In the past few years, several groups have used the FP to investigate the evolution of the 
mass-to-light ratio of held spheroidal galaxies either using direct measurements of velocity 
dispersion (Treu et ah 1999, 2001a,b, 2002, 2005; van Dokkum et ah 2001; van Dokkum & 
Ellis 2003; Bernardi et al. 2003; Gebhardt et al. 2003; van der Wei et ah 2004) or estimating 
stellar velocity dispersions from the separation of multiple images of gravitationally lensed 
background objects (Kochanek et al. 2000; Rusin et al. 2003; van de Ven et al. 2003; Rusin 
& Kochanek 2005). The Lenses Strncture and Dynamics (LSD) Survey (Koopmans & Treu 
2002,2003; Treu & Koopmans 2002,2003,2004; hereafter collectively KT) measured stellar 
velocity dispersions of gravitational lenses thus bridging the gap between the two methods. 
For reference, previous measurements are snmmarized in Table 2. 

The hrst thing to notice from the table is that all published measurements indicate faster 
FP evolntion than that fonnd for cluster galaxies (—0.49 ± 0.05, van Dokknm et al. 1998) 
although in some cases the diherence is not statistically signihcant. The second important 
point is that most measurements, even those based on lensing properties, agree with one 
another to within the quoted 1-a uncertainties. This can be explained by the strnctural 
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homogeneity of early-type galaxies which appear to have close to isothermal mass density 
prohles (e.g. Rusin et al. 2003b) as a result of a still unexplained dark- luminous matter 
conspiracy (e.g. Treu & Koopmans 2004). 

If we take the measurement from this paper as the reference, all measurements listed 
in Table 2 are in agreement within the uncertainties except for van Dokkum & Ellis (2003) 
and Rusin et al. (2003) who measure a slightly slower evolution, and Gebhardt et al. (2003), 
whose work appears to indicate faster evolution. The latter discrepancy might be explained 
by a selection bias of the type shown in Figure 15, particularly since the Gebhardt et al 
study used considerably shorter exposure times at comparable redshifts. Other differences 
might be explained by the different mass ranges in the various samples, contamination by 
bulge-dominated spirals in shallower imaging data, or simply by small number statistics and 
different htting techniques. If mass-to-light ratios follow the tracks illustrated in Figure 21, 
as a result of a series of secondary bursts overlayed on an old stellar population, the aver¬ 
age evolutionary rate will depend on the fraction of objects caught in the active luminous 
phase. Similarly, different htting techniques could be more or less sensitive to “outliers” and 
yield varied answers and error estimates. In short, small samples will be very sensitive to 
huctuations in this number and to the adopted htting technique. 

Previous studies have often emphasized the diherences between the results, arguing 
whether or not held early-type galaxies have younger stellar populations than their counter¬ 
parts in clusters. Our sample allows us to hrmly establish that held early-type galaxies have, 
on average, younger stellar populations than cluster galaxies, as a result of a more diverse 
star formation history, i.e. a wider range in epochs of star formation and secondary activity 
extended at relatively recent cosmic times. 


5.6. Independent diagnostics of recent star formation 

In T05 we showed that the mass-to-light ratio measured from the FP correlates well with 
the rest frame B-V color, indicating that deviations from the passive trend (dehned in that 
paper by that observed for cluster early-types; and measured with the parameter SAM/Lb, 
Section 5.1) is associated with younger stellar populations. In this Section, we will extend 
this analysis, using spectral ([O II] and 11(5) and morphological (blue cores) diagnostics as 
independent probes of the star formation history of our sample of E-I-SO galaxies. We will 
use [O II] as an indicator of ongoing star formation (e.g. Kennicutt 1992), since the life time 
of the hot OB stars required to produce photo-ionizing photons is very short. Similarly, as 
strong H(5 absorption occurs primarily in A stars with main sequence lifetimes of 1 Gyr, 
we use its occurrence to measure the recent activity. 
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Prior to using blue cores and emission lines as star formation diagnostics, we must 
consider what fraction of this activity might be due to non-thermal AGN. For this purpose 
we identify nuclear activity on the basis of X-ray luminosity as discussed in Section 4. 
Adopting a conservative threshold of Lx > 10^^ erg s“^, we hnd only two of the 14 galaxies 
with blue cores are likely to be AGN dominant, with a further a possible LLAGN. Similarly, 
three strong [O II] emitters (i.e. EW> IsA; 2 early-types, 1 Sa-l-b) without blue cores are 
identihed as AGNs. These objects are excluded from the analysis in the rest of this section, 
guaranteeing that the remaining blue cores and [O II] are indeed due to star formation. 

To visualize the connection between star formation indicators, morphology and mass to 
light ratio, we plot in Figure 22 the distribution of the offset from the redshift-dependence of 
the FP evolution of massive cluster spheroidals SAM/L^ (Equation 11). Galaxies with signs 
of recent star formation (blue-cores, strong [O II] emitters, and H(5 equivalent widths< —SA) 
show in general a larger offset than the quiescent population. However, even the quiescent 
population evolves more rapidly than the cluster population {SAM/Lb = 0). Quantitatively, 
for the E-I-SO galaxies, the average SAM/ Lb is —0.23±0.02 for the entire sample, —0.39±0.07 
for the strong galaxies, —0.39 ± 0.10 for the star forming blue cores, and —0.17 ± 0.12 
for the star forming [O II] emitters. 

In conclusion, galaxies whose images and spectra reveal evidence for recent star forma¬ 
tion have lower mass to light ratios at a given redshift than their quiescent counterparts. 
No signihcant difference is found for the sample of galaxies with [O II] emission, perhaps 
because of small number statistics (11 E-|-S0s, excluding the 2 AGN) or perhaps because of 
the different timescales ([O II] is a measure of instantaneous star formation, while M/Lb 
is a measure of integrated star formation). However, the fraction of [O II] E/SO emitters 
(11/54) for which the line is visible with our spectroscopic setup is signihcantly higher than 
for early-type galaxies in the local field (~ 0.05 Treu et ah 2002), providing further evidence 
for enhanced star formation. 

Studying the evolution of H(5 absorption we can obtain an independent check on our 
inferences based on the evolution of the mass to light ratio and perform an independent test 
of the secondary activity postulated in Section 5.3.1. In Figure 23 we show H5 equivalent 
width vs. redshift for all E-I-SO galaxies ^ and compare it with the predictions of our simple 
(thick red line) and composite population models (black dashed lines). We hnd that the 
distribution of rest-frame H(5 equivalent widths is inconsistent with any single burst model 
formed at high redshift and recent star formation activity is required. As for the M/Lb, the 


^Note that H5 is only visible for a subset of objects (74/163 E-|-S0s), i.e. when redshifted into the 
wavelength range covered by our spectrograph. 
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distribution is well-reproduced by a composite population where the bulk of stellar mass is 
formed at high redshift (e.g. Zfi ~ 3) and on average 10 % of the stellar mass is added at a 
later time {zf 2 ~ 0.5 — 1) in secondary bursts. 

Evidence for downsizing in the population, discussed initially by T05, is also seen in 
the H(5 absorption. To illustrate this, coadded spectra of all the E-I-SO galaxies binned in 
redshift and mass are shown in Figure 24. For galaxies with masses smaller than 10^^ Mq 
we hnd an average H(5 equivalent width ~ —3 A, inconsistent with the expected value for an 
old stellar population (c.f. Figure 24; see also van Dokkum & Ellis 2003 calculate -lA 
for stellar populations formed at 2 : ~ 3 and observed at z ^ 0.9). 


6. Summary and conclusions 

We have obtained ultra-deep spectroscopy of 163 held E-I-SO and 61 bulge dominated 
Sa-|-b galaxies, selected morphologically from the GOODS-N held to a magnitude limit of 
zg < 22.43. Redshifts are measured for all objects in the sample, while stellar velocity 
dispersions are obtained for 141/163 E-|-S0s and 40/61 Sa-|-bs. The E-|-S0s with measured 
stellar velocity dispersion represent an unbiased and representative subsample of the parent 
morphological population. Structural parameters are obtained by htting prohles to the 
ACS images. Morphological asymmetries and blue cores are identihed on the basis of the 
multicolor ACS images. Blue core spheroidals represent approximately 10% of the population 
(14/163). When visible with our spectroscopic setup, [O II] and H(5 are measured and used 
as secondary star formation indicators. 

Our redshift catalog is cross-correlated with the deep 2Ms X-ray catalog of the Chandra 
Deep Field North (Barger et al. 2003) to identify AGN and star forming galaxies. Nine 
E-|-S0s, including 2 blue cores, and 6 Sa-|-bs are identihed as AGN based on an X-ray 
luminosity above 10^^ erg s“^. 21 additional galaxies (16 E-|-S0s, 4 Sa-|-bs, IS) are also 
detected in the X-ray with luminosities between 10^° and 10^^ erg s“^, which are interpreted 
as mostly due to a low luminosity AGN or to ongoing star formation. The vast majority 
(11/14) of blue core E-|-S0s are not detected in the deep Chandra exposures, indicating that 
the blue cores are the result of spatially concentrated star formation activity, and not of non 
thermal AGN emission. 

We study the evolution of the Fundamental Plane as a diagnostic of the star formation 
and mass assembly history of early-type galaxies, using Montecarlo simulations to rigorously 
control selection ehects. The main results can be summarized as follows: 

1. The mass-to-light ratio of held spheroidals evolves on average as d\og{M/L b)/ dz = 
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—0.721q;q 5±0.04, i.e. faster then that of cluster spheroidals (—0.49 ±0.05 van Dokkum 
et ah 1998). This is consistent with the stellar population of held spheroidals being 
on average younger. Although the difference between the FP evolution of cluster and 
held spheroidals (Table 2) is in the sense expected by hierarchical models (Diaferio et 
ah 2001; De Lucia et al. 2004), quantitatively it is not as pronounced as expected (van 
Dokkum et al. 2001; Treu 2004). 

2. The evolution of the Fundamental Plane, and thus of the ehective mass to light ratio, 
depends strongly on mass: the less massive the galaxy, the faster the evolution. The 
mass-to-light ratio of the most massive galaxies (> IO^^'^Mq) evolve as that of an old 
stellar population {zf > 2), with an upper limit of ~ 1% to the fraction of stellar 
mass that can have formed at z < 1.2. The evolution of the mass to light ratio of the 
less massive galaxies is inconsistent with that of a single burst stellar population. A 
signihcant fraction of the mass (up to 20-40 % for masses below 10^^ Mq) must have 
formed more recently than Zf 2 = 1.2 consistent with signihcant recent growth of this 
population. 

3. The intrinsic scatter of the FP is found to increase with redshift, after allowing for 
evolution of the slopes, as dims (log M / L^) / dz = 0.040 ± 0.015, consistent with the 
idea that we have approached an epoch of starformation activity in this population. 

4. The scenario of mass dependent star formation history and signihcant recent growth for 
E±S0 with mass below 10^^ Mq is supported by independent star formation indicators. 
The mass-to-light ratio of spheroidal galaxies with blue rest frame colors (discussed 
in a companion paper, T05), strong emission or balmer absorption ([O II] and H5), 
and morphological anomalies (blue cores) is found to evolve more rapidly than that of 
quiescent galaxies. The fraction of E±S0s galaxies with [O II] emission is larger (13/54, 
or 11/54 excluding two emitters identihed as AGN based on the X-ray emission) than 
in the local universe (~0.05; Treu et al. 2002 ). The distribution of Hh absorption for 
our sample is well reproduced by the same two-bursts models describing the evolution 
of the FP. 

We promote, therefore, a picture in which all spheroidals except the most massive ones 
continue to assemble and form young stars via subsequent episodes of star formation. The 
associated rejuvenation is sufficient to signihcantly inhuence the location of the galaxy on the 
FP and, of course, its colors. As a result, a fraction of spheroidal galaxies will have blue colors, 
and be missed in color selected samples (e.g. Jimenez et al. 1999). Likewise, passive models 
will be an inadequate representation of the evolutionary behavior. Introducing secondary 
bursts implies that for a given luminosity function (or equivalently observed surface density; 
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Im et al. 2002; Cross et al. 2004) will lead to an increased evolutionary rate in the number 
density evolution (Treu 2004). 

Ultimately we seek to connect the evolution of the mass to light ratio for various masses 
with the broader issue of how the number of spheroidals is growing with cosmic time (c.f. 
Franx 1993). To illustrate a possible approach to this important larger picture, we explore 
combining the FP measurements made in this paper with constraints on the evolution of 
the luminosity function (LF), for example morphologically-selected E-|-S0s from the WFPC2 
survey discussed by Im et al (2002). 

Im et al. (2002) chose to represent the evolution of their E-I-SO LF with two parameters: 
Q describing the luminosity evolution (equivalent to —2.bd\og{M/L-Q)/dz), and m describing 
the number density evolution, i.e. the characteristic density parameter in the Schechter 
(1976) luminosity function which is postulated to evolve as {1+z)'^. With luminosity function 
data alone, there is a degeneracy betweem Q and m: stronger luminosity evolution implies 
a stronger number density evolution. Using our measured value for the luminosity evolution 
inferred independently from the FP (equivalent to Q = 1-80 ±0.15) in combination with the 
surface shown in their Figure 24, we eliminate the degeneracy and obtain m = —0.6 ±0.5. 

This value of m implies an abundance of morphologically selected E±S0s at ~ 1.2 of 
60^20 ^ of local value. If correct, the decline in the number density of red galaxies (Bell 
et al. 2004) would be explained in part by an evolution of the number density of E±S0s 
galaxies (at least those morphologically identihable, e.g., van Dokkum & Franx 2001) and 
in part by their blue colors. Such an interpretation would apply equally to samples of color- 
and morphologically- selected higher redshift E±S0 galaxies (Menanteau et al. 1999; Treu 
& Stiavelli 1999; Stanford et al. 2004). 

Although a crude resolution of the differences seen for the various datasets, largely be¬ 
cause of the uncertainties in combining different datasets and the limited size (and thus 
cosmic variance noise) of the Im et al. 2002 sample, by combining of our precise FP mea¬ 
surements with more extensive luminosity functions a more robust empirical picture of the 
evolution of E±S0s in the last 9 Gyrs can be determined. 
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Table 1. Catalog of spheroidals and bulge dominated galaxies 


ID 

RA 

DEC 

T 

t>4 re4 

V6 

^■ee 

is 

»’e8 

Z9 

»’e9 

z 


[Oil] 

H5 

S/N 

texp 

1287 

189.13248 

62.15182 

0 

_ 

21.90 

0.81 

20.80 

0.45 

19.95 

0.47 

0.8457 

319±22 

-0.4±0.3 

1.3±0.2 

22 

19800 

648 

189.16547 

62.13849 

0 

22.30 0.76 

21.21 

0.37 

20.58 

0.29 

20.24 

0.28 

0.2469 

65±13 

- 

- 

26 

19130 

609 

189.19130 

62.10799 

0 

- 

22.45 

0.74 

21.09 

0.55 

20.51 

0.49 

0.7976 

247±16 

5.4±0.6 

-2.2±0.4 

12 

19130 

1286 

189.13710 

62.15260 

2 

- 

22.41 

1.21 

21.06 

0.85 

20.36 

0.76 

0.8461 

231±16 

-0.6±0.6 

1.7±0.4 

12 

19800 

655 

189.19568 

62.13973 

1 

- 

21.91 

0.78 

20.72 

0.81 

20.35 

0.75 

0.6421 

99±7 

- 

3.7±0.4 

14 

19130 

1253 

189.13830 

62.14290 

3 

- 

22.23 

0.03 

21.90 

0.03 

21.24 

0.10 

0.9349 

- 

3.8±0.2 

1.8±0.2 

18 

19800 

1236 

189.13470 

62.13730 

0 

- 

22.91 

0.61 

21.59 

0.40 

20.92 

0.35 

0.8503 

202±11 

0.5±0.4 

-2.0±0.3 

14 

37800 

1278 

189.15150 

62.14850 

2 

- 

22.51 

0.14 

21.57 

0.12 

21.19 

0.12 

0.5125 

108±7 

- 

- 

22 

37800 

681 

189.16202 

62.15102 

1 

- 

23.28 

0.40 

21.87 

0.34 

21.17 

0.32 

0.8422 

318±28 

-0.7±0.5 

-0.3±0.3 

13 

38930 

635 

189.20297 

62.12884 

2 

- 

23.24 

0.59 

21.73 

0.45 

21.05 

0.41 

0.8195 

187±16 

-0.5±0.5 

-2.1±0.4 

11 

38930 


Note. — For each object we list ID, RA, DEC, morphological type, AB magnitudes (from fitting) and effective radii 

(in arcsec), redshifts, measured velocity dispersion crap in km s~^, [O II] and equivalent widths (in A, rest frame), average 
signal-to-noise ratio (A~^; observer frame) and spectroscopic exposure time (in seconds). Photometric parameters through the 
64 filter are derived only for objects at 2: < 0.34. Formal galfit uncertainties on magnitudes and effective radii are 0.01-0.02 mags 
and a 1-10% in effective radii. Total uncertainties on the combination of photometric parameters that enters the Fundamental 
Plane log Re — 0.32SBe are 0.02-0.03. The Table is published in its entirety online, the first 10 entries are shown here for 
guidance. 



Table 2. 


Summary of previous field measurements 


dlog{M / Lb) / dz 

buedshift 

sample 

Source 

Notes 

-0.50±0.19 

O.O-l.O 

22 

Rusin &: Kochanek 2005 

lensing-general; lenses 

-0.54±0.06 

0.6-1.0 

9 

van Dokkum & Ellis 2003 

direct; morphology and magnitude 

-0.56±0.04 

O.O-l.O 

27 

Rusin et al. 2003 

lensing-isothermal; lenses 

-0.59±0.15 

0.2-0.6 

18 

van Dokkum et al. 2001 

direct; morphology and magnitude 

-0.62±0.13 

O.O-l.O 

26 

van de ven et al. 2003 

lensing-isothermal; lenses 

-0.71±0.19 

0.6-1.1 

6 

van der Wei et al. 2004 

direct; morphology, magnitude and color 

-0 72+^-^® 

'^-0.11 

0.1-0.7 

30 

Treu et al. 2002 

direct; morphology, magnitude and color 

-0.72±0.10 

0.5-1.0 

5 

Treu & Koopmans 2004 

direct; lenses, early-types only 

-0.75±0.17 

0.5-1.0 

5 

Treu & Koopmans 2004 

lensing-isothermal; lenses, early-types only 

-0.80..-96 

0.3-1.0 

36 

Gebhardt et al. 2003 

direct; magnitude and absorption lines 

-0.72t°;°^ ± 0.04 

0.2-1.2 

141 

this work 

direct; morphology and magnitude 


Note. — In the notes column, “direct” refers to spectroscopically determined velocity disper¬ 
sions, “lensing-isothermal” refers to velocity dispersion estimate from image separation using an 
isothermal model, “lensing-general” refers to velocity dispersion estimate from image separation 
using a more general model (see Rusin &: Kochanek 2005); the measurements for Gebhardt et ah 
2003 are the average brightening to z = 0.8 and z = 1.0, respectively. 
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A. Transformations from ACS magnitude to rest frame B and V magnitudes 

Throughout this paper, rest frame Landolt B and V magnitudes are obtained from the 
two ACS hlters which are closest to the redshifted B and V central wavelengths (i.e. 4405 and 
5470 A, e.g., Fukugita et ah 1996). Given the redshift interval considered here 0 < z < 1.25, 
a pair of ACS hlters is always sufficient to obtain rest frame B and V magnitudes within few 
hundreds of magnitudes. Hence, the transformations adopted here can be expressed in the 
form: 


B = Ue + aBvbv(&4 - ^' 6 ) +/^Bvbv - dm [0 <= < 0.34] (Al) 

H = is + aBivi (^'6 - ^s) +/^Bivi - dm [0.34 < 2 : <= 0.76] (A2) 

B = Zg + aBziz(i 8 - Zg) + /^Bziz “ DM [0.76 <= < 1.25] (A3) 

H = ue + Q;vvbv(&4 - ^' 6 ) +/Svvbv - dm [0 <= 2 : < 0.07] (A4) 

H = is + avivi (^'6 - ^s) +/Svivi - dm [0.07 < 2 : <= 0.41] (A5) 

H = 2:9 + avziz(i 8 - 2 : 9 ) +/^Vziz - dm [0.41 <= 2 : < 1.25] (A 6 ) 


where the coefficients a and (3 are a function of redshift, while the bolometric distance 
modulus DM is a function of redshift and cosmography. Polynomial expansions of the 
coefficients are given in Table 3. 
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N 

O^Bvbv 

/^Bvbv 

O^Bivi 

/^Bivi 

O^Bziz 

P'Bziz 

0 

0.824411 

-0.00169895 

2.44316 

0.748236 

1491.01 

-796.805 

1 

-3.29360 

1.22223 

-4.42301 

-3.72483 

-8947.85 

4776.19 

2 

3.40645 

-8.36258 

1.53361 

7.43317 

22264.2 

-11813.0 

3 


18.9313 


-3.58153 

-29340.5 

15426.2 

4 





21564.4 

-11200.7 

5 





-8367.60 

4278.26 

6 





1336.45 

-669.418 

N 

Q^Vvbv 

P'Vvhv 

Q^Vivi 


^Vziz 

/^Vziz 

0 

0.186741 

-0.0274962 

1.42588 

-0.0394309 

-77.0651 

46.5707 

1 

-2.62562 

0.995802 

-5.13394 

1.23626 

668.581 

-399.767 

2 

7.19822 

0.0925006 

4.22521 

-1.42416 

-2294.21 

1398.85 

3 




1.88833 

4056.19 

-2539.31 

4 





-3936.67 

2531.96 

5 





1999.86 

-1319.01 

6 





-417.432 

281.589 


Table 3: Taylor expansion coefficients of a and (3 aronnd redsliift z = 0. Redsliift intervals 
of validity are given in Eqnations A1 to A6 



36 


B. X-ray-based classification scheme 

In this paper we adopted a simplified classification scheme, based on X-ray Inminos- 
ity, to identify galaxies with nuclear activity (XAGN), or with low luminosity nuclear ac¬ 
tivity and/or ongoing star formation (XSF/LLAGN). The working dehnitions of AGN or 
XSF/LLAGN - admittedly schematic - are motivated by local studies of early-type galaxies 
which show that typically an AGN is required to power a luminosity above 10^^ erg s“^ 
(with the exception of cluster brightest cluster galaxies - not present in our sample - where 
the IGM can reach these luminosities e.g. Pellegrini & Giotti 1998, O’Sullivan, Forbes & 
Ponman 2001). At fluxes below 10^^ erg s“^ the main contributors to the X-ray luminosities 
are LLAGNs, star formation related phenomena such as HMXB, low-mass X-ray binaries 
(LMXB), and emission from the diffuse plasma, with tipical temperatures below 1 keV. 

We can estimate the contribution of the two latter mechanisms as follows. The con¬ 
tributions from LMXB is proportional to the stellar mass and thus scales linearly with the 
B-band luminosity and it is generally small with respect to the total X-ray emission ob¬ 
served in local early-type galaxies (see O’Sullivan, Forbes & Ponman 2001; Kim & Fabbiano 
2004). Using the relation from Kim & Fabbiano Lx (LMXB) = (0.9±0.5) 10^° erg s“^ Lb / 
Lb,o and our measured Lb we hnd that for our sample^ the median contribution of LMXBs 
to the total X-ray emission is 16±9%. 

As far as diffuse plasma emission is concerned, we can conservatively estimate its con¬ 
tribution by using the empirical scaling of the total soft X-ray luminosity with Lb as mea¬ 
sured for E-I-SO galaxies in the local universe, excluding AGN hosts and brightest cluster 
galaxies. Since local early-type galaxies are mostly quiescent, the excess with respect to 
the local relationship provides an estimate of the contribution due to star formation or 
LLAGN in high redshift sources. For example O’Sullivan, Ponman & Gollins (2003) hnd 
Lo. 5-2 — 4 X 10^^ X (Lb/10^^)^'’^, albeit with a scatter of an order of magnitude. Applied 
to our sample, this relationship results in an expected X-ray luminosity corresponding to 
6±4% of the total observed one. As an example of the uncertainty in this estimate, adopting 
instead the relationship in Eq.l from O’Sullivan, Forbes & Ponman (2001) raises the fraction 
to 25 ± 15%. We conclude that in general the X-ray luminosity of our distant E-I-SO with lu¬ 
minosities in the range 10^°-10^^ erg s“^ is due largely to star formation and/or LLAGN and 
cannot be explained entirely by the plasma and LMXB emission expected for local E-I-SO. 
However, given the high intrisic scatter of the local relationships, XSF/LLAGN might not 
be needed in some individual cases. 


•^logLe/l^B,© is in the range 9.2-10.8 average 10.2, after evolving to 2 : = 0 using the evolution of the FP 
measured in Section 5. 
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Fig. 1.— Comparison of morphological classifications based on the vO.5 and vl.O releases of 
the GOODS data, referred to as “old” and “new” respectively. Types T are assigned based 
on the following scheme: T=-2=star, -l=compact, 0=E, 1=E/S0, 2=S0, 3=Sa+b, 4=S, 
5=Sc+d, 6=Irr, 7=Unclass, 8=Merger, 9=Eault (Abraham et ah 1996; Treu et ah 2003). 
Left panel refers to all classes, right to spheroidals only. 


Fig. 2.— Color-composite images (R= 2 : 9 , 0=^8, B=n6) of all early-type galaxies in the Keck 
spectroscopic sample, sorted by redshift. The color scheme is chosen so that a local early- 
type galaxy would appear white. Images are 3" on a side. For each galaxy redshift, ID, and 
morphological type are shown in the upper left, lower left and lower right corner respectively. 
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Fig. 2.— Color-composite images of all early-type galaxies (continued). 


Fig. 2.— Color-composite images of all early-type galaxies (continued). 


Fig. 3.— Color-composite images (R= 2 ; 9 , 0=^8, B=n6) of all the early-spirals in the Keck 
spectroscopic sample, sorted by redshift. 
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Fig. 4.— Examples of de Vaucouleurs fits: for galaxy types —1 < T <4 (one per row), for 
three different magnitude bins, ~ 20, ~ 21, ~ 22. For each galaxy the original image 

is shown alongside residuals from the best £t model. Notice the disk-like residuals in 
the SO galaxies and the spiral arm residuals in the Sa-fb galaxies conhrming the vl.O visual 
classification. Redshift, magnitude and ID are listed for each object. Each panel is 3'.'0 x 2'.'7. 


Fig. 5.— Montage of ’blue core’ early type galaxies, selected as defined in the text, sorted 
by redshift. Each image is 3” on a side. Only the two highest redshift examples appear to 
host active galactic nuclei, as indicated by the large X-ray luminosity (> 10"^^ erg s“^). The 
bulk of the blue light arises from recent star formation. 
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Fig. 6.— Distribution of colors (panels a and b) and zg magnitudes (panel c) for the full 
sample of z<^ <22.43 E+SO galaxies classified in the hnal vl.O-zcatalog (empty histogram) 
compared to the E+SO sample (vl.O morphologies) targeted for spectroscopy (solid his¬ 
togram). The sampling rate is approximately independent of color, but drops off beyond 
z^ ~ 21.5 as discussed in Section 3.1. 
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z-0.465 z,-20.01 ID-244 z-0.512 z,-21.19 ID-1278 z-0.475 z,-21.72 ID-243 



z=0.766 z,=20.27 ID=774 z=0.850 z,=20.92 ID=1236 z=O.BOl z,=21.72 ID=1076 



2=1.013 z,=20.29 ID=761 2=1.140 z,=20.62 1D=1373 2=0.960 2,=21.B1 ID=1761 



Fig. 7.— Examples of flux-calibrated spectra - smoothed with a boxcar hlter of width 
gA- for three different magnitude bins, Zg ~ 20, Zg ~ 21, zg ~ 22 and redshifts {z ~ 
0.5,2: ~ 0.8,2: ~ 1). The signal-to-noise ratio decreases with magnitude and exposure time 
as expected. Fluxes are in units proportional to photons per second. 
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Fig. 8.— Redshift distribution of the early type spiral and E+SO sample. Spectra with 
S/N>8 A-i in the observers frame, for which velocity dispersions are typically available, are 
shown as a hlled histogram. 
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Fig. 9.— Example of an intermediate S/N spectrum (top) and the highest redshift spectrum 
with measured velocity dispersion (bottom). Both spectra are boxcar smoothed with a bA 
hlter. The red line is the best htting stellar template smeared at the best htting velocity 
dispersion. Hatched regions represent those masked out during the £t. The top scale shows 
the rest frame wavelength and the derived quantities. 
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Fig. 10.— Success rate of velocity dispersion measurements. The upper panel shows the 
redshift distribution of early-type galaxies in the spectroscopic sample (empty histogram) 
and those with measured velocity dispersions (solid histogram). The middle panel shows the 
equivalent Zg luminosity distribution of the two samples. The lower panel shows the <S/N> 
(A“^;observer frame) distribution of the two samples (note the inverted logarithmic x-axis). 
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Fig. 11.— Velocity dispersion tests: I. Solid points represent comparisons with LRIS mea¬ 
surements from van Dokkum & Ellis (2003), open points represent internal DEIMOS com¬ 
parisons for objects observed on different mask and slitlets configurations. 
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Fig. 13.— Location of survey galaxies in FP-space (the projection corresponds approximately 
to an-edge on view of the local FP) binned in redshift. Solid pentagons represent E+SO 
galaxies with a >= 100 km s“^, open pentagons represent E+SO galaxies with a < 100 
km s“^, while starred squares represent Sa+b galaxies. The local relationship is shown as a 

solid line. 
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Fig. 14.— Offset from the local FP for the early-type and early-spiral galaxies, coded by 
mass (bigger symbols represent bigger dynamical masses). Notice that the offset from the 
local relationship and the scatter decrease with mass. 
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Fig. 15.— Effects of magnitude selection on the determination of the FP. The panels and 
symbols are as in Figure 13 but the FP is now projected on the M-M/L plane. The local 
relation is shown as a dotted line for comparison. Hatched regions are excluded by our 
Zg < 22.43 magnitude limit. The more densely hatched regions indicate the selection limit 
applicable for the lower end of the redshift bin, the more lightly hatched one to the upper 
end (assumed at z=l.l for panel e). The plots are for illustrative purposes only; the limits 
are approximate to a few hundreds of a magnitude because a hxed /c-color correction was 
adopted and because the model photometry used to construct M/Lb is somewhat different 
from the SExtractor photometry in the selection process. The true selection criteria in the 
Zg band are faithfully implemented in our Montecarlo-Bayesan method. 
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Fig. 16.— Probability contours for the evolution of the intercept d'y/dz and maximal evo¬ 
lution of the scatter da.~f/dz (dashed=68%; solid=95%). The maximum of the probability 
distribution is indicated by the solid pentagon. Thin contours represent the same probabil¬ 
ity levels for a smaller sample of 119 E-|-SOs, where 21 SOs have been excluded as possible 
early-spirals contaminants, see Section 2.2. The corresponding values in terms of M/L are 
shown on the top and right axis. 
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Fig. 18.— Bottom panel: the shaded area represents the 68% probability contonr for the 
evolntion of M/L as a fnnction of M, obtained from Fignre 17 with Eqnation 9. The solid 
red line corresponds to the best htting model. Top panel; distribntion of dynamical masses 
for onr sample. 
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Fig. 20.— Redshift of formation as a function of evolution of the effective mass to light ratio 
for a set of simple stellar population models. 
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Fig. 21.— The solid black lines show the evolution of the mass-to-light ratio for a family of 
two-bursts models, where the majority of the stellar mass is formed at Zfi = 5 and the rest 
(1-5-10% as indicated on the plots) is formed at Zf 2 < Zfi (secondary bursts are shown in 
steps of dz = 0.1, in the range z = 0.3 — 3). For illustration the evolution of the mass-to-light 
ratio of the old component is also shown as a thick solid red line. The observed offset for our 
sample of E-I-SO galaxies (coded by mass as in previous figures), is also shown for comparison. 
Note that this simple model covers the observed range of parameters, demonstrating that a 
simple two burst model is sufficient to reproduce the observations. 
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Fig. 23.— H5 as a function of redshift for a family of two bursts models, consisting of a 
main burst forming 90% of the stellar mass at Zfi = 3, plus a secondary burst forming 10% 
of the stellar mass at Zf 2 < Zfi (lines as in Figure 21; models from Bruzual & Chariot 2003; 
solar metallicity). Measurements of H(5 for our sample of E+SOs galaxies are shown as hlled 
pentagons. 
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Fig. 24.— Co-added spectra of E-I-SO galaxies, binned in mass and redshift. H(5 in absorption 
is seen to be stronger for smaller mass systems, consistent with younger stellar ages than 
their massive counterparts. 




